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ABSTRACT 

We show here a component of the meridional circulation developing at medium-high latitudes (40- 
50°) before the new solar cycle starts. Like the torsional oscillation of the zonal flows, this extra 
circulation seems to precede the onset of magnetic activity at the solar surface and move slowly 
towards lower latitudes. However, the behavior of this component differs from that of the torsional 
oscillation regarding location and convergence towards the equator at the end of the cycle. The 
observation of this component before the magnetic regions appear at the solar surface has only been 
possible due to the prolonged solar minimum. The results could settle the discussion as to whether 
the extra component of the meridional circulation around the activity belts, which has been known 
for some time, is or is not an effect of material motions around the active regions. 
Subject headings: Sun:helioseismology, Suniinterior, Sun:meridional circulation, Sun:activity 



1. INTRODUCTION 

Meridional circulation has became an important 
layer in the fiux- transport solar-dynamo models 
Miesch I 12005. and refe rences therein). Until the de- 
velopment of local helioseismology methods, the ob- 
servation of meridional circulation was limited to 
the surface layers by tracing^ mag;netic elements 
(lKomni._H oward fc Harvev 1993; Snodgrass fc Dailev I 

119961 : iNesme -Ribes, Meunier & Vince 199't) or by 

using surfa c e Do ppler velocities (|Hathaw^n Il996f ). 
iGiles et al. I (|l997l ) showed the first meridional circula- 
tion inferences under the solar surface by applying a par- 
ticular local helioseismology technique, time distance. 
Although the full picture of the meridional circulation 
throughout the convection zone remains elusive, and re- 
cent studies suggest that very long se ries of data are 
required to infer the flows deep down (|Braun fc Birch I 
[2OO8), local helioseismology has been able to give de- 
tailed information on these flows in the subsurface lay- 
ers. In particular, the variation wit h the solar cycle has 
been the focus of several studies ("Haber et al. ' '20021; 
iBasu fc Antial l2003l : Gonzalez Hernandez et al. 200i 
IZaatri et al. I I2OO6I : [Gonzalez Hernandez et. all i2008t ). 
which have shown that the overall amplitude of the 
meridional circulation is anticorrelated with the mag- 
netic activity. This trend has been confirmed by 
tracing weak magnetic eleme nts on the solar surface 
dHathawav fc Rightmire]l2Q10|). 

Chou and Dai (200T|) presented the first evidence of 
an additional subsurface component of the meridional 
circulation that varied with the solar cycle located 
around the center of the sunspot distribution. Apply- 
ing the time-distance technique to observations from the 
Taiwan Oscillation Network (TON), they identified a 
residual of the meridional flow that correlated with the 
location of the surface activity for the declining phase 
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of s olar cycle 22 and the increasi ng phase of solar cycle 
23. iBeck. Gizon fc Duvall I (|2002f ) confirmed the discov- 
ery using Michelson Doppler Imager (MDI) data from 
1996 to 2001, and suggested this extra component could 
be a new component of the solar dynamo. In both cases, 
the authors conclude that this new component could be 
related to the well known phenomena of material moving 
inward/outward from active regions, although they no- 
ticed that the inferen c es sug gested a larger than expected 
depth range. Gizon ' ('2003') tested the possibility of the 
extra circulation being related to the inflows around large 
active regions recorded by local hel i oseismology analysis 
(I Zhao. Kosovichev fc Duvah I I2OOII: iHaber et aTl 120041 : 



Kom m et al. 1120041 : iBraun. Birch fc Lindsev Il2004f ). He 

analyzed two Carrington rotations eliminating the con- 
tribution from areas with surface magnetic activity and 
found that the extra component decreased. 

The continuous operations of the Global Oscillation 
Network Group (GONG) throughout the full solar cycle 
23 has given us the unprecedented opportunity to study 
the subsurface meridional circulation without interrup- 
tion. Using these data, Gonzal ez Hernandez et. all 
(|2008l ) revealed that even after agressively masking sur- 
face magnetic activity over a 5.5-year period (2001-2006), 
the extra component of the meridional circulation re- 
mained present, although attenuated, under the solar 
surface. Here we extend this analysis to the prolonged 
solar minimum and find that the formation of the bumps 
at medium latitudes precedes the onset of the activity, 
confirming the nature of this phenomenon as indepen- 
dent of the surface manifestation of active regions. 

2. DATA ANALYSIS 

We extend the s tudy described by 
iGonzalez Hernandez et. aTl (|2008f ) to years 2007, 
2008 and 2009 to include the recent prolonged solar 
minimu m. In s ummary, we apply standard ring-diagram 
analysis (jHillJ'l988), a local helioseismology technique, 
to GONG high-resolution Dopplergrams to infer the 
meridional flow from the solar surface to a depth of 
approximately 16 Mm. 
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Fig. 1. — Yearly averages of the meridional flow obtained by ring-diagram analysis of continuous GONG data at four different depths. 
The variation with the solar cycle clearly observed at the superflcial layers is less pronounced at deeper layers. The extra circulation 
(bumps) is also clearly visible in the shallow layers. 



The ring diagrams method uses medium/high-degree 
waves propagating in locahzed areas of the Sun to ob- 
tain an averaged horizontal- velocity vector for that par- 
ticular region. The input data used for the analysis are 
high resolution Dopplergrams. In the typical analysis, 
a particular area of the Sun (16° square, apodized to 
15°-diameter) is tracked and remapped, and the analy- 
sis of the corresponding tridimensional power spectrum 
of solar oscillations renders an average horizontal v eloc- 
ity ve ctor for the area at different depths (Corbard e t al. I 
l2QQ3f ). By analyzing a mosaic of these patches, it is possi- 
ble to infer a three-dimensional velocity map in the depth 
range where the waves propagate. Typical ring-diagram 
analysis uses 1664-minute series of data with a resolution 
of about 1.5 Mm per pixel at the center of the disk. We 
obtain horizontal velocity maps from disk center to ap- 
proximatel y 52.5° in each direc tion with patches centered 
every 7.5° (|Haber et al. |[2QQ2l ). To study the meridional 
circulation we isolate the Vy component of the calculated 
flows and use patches centered at no more than 22.5° 
heliocentric longitude to minimize projection effects. 

Figure [1] presents the yearly- averaged meridional flows 
calculated for the period of August 2001 to April 2009. 
The results show the well known variation of the over- 
all amplitude of the meridional circulation with the solar 
cycle: smaller amplitudes in 2001-2002 (maximum activ- 
ity) and larger amplitudes towards 2008-2009 (minimum 
activity). The variation with the solar cycle is less pro- 



nounced at deeper depths, but it i s still visible at high 
latitudes. Gizon fc Rempel I (|2008f ) present a model to 
account for independent observations of the meridional 
flow at the surface and at a depth of 60 Mm. In such 
model, the meridional flow at these two depths are an- 
ticorrelated. Although our analysis only returns infor- 
mation over the upper 15 Mm, the depth dependence of 
the results is such that it would be consistent with their 
findings, with a reversal in the time varying component 
at depths below our accessible range. 

Fitted low order polynomials to 60-day averages of 
the meridional flow at a depth of 5.8 Mm are shown in 
Figure [2l Positive velocities are taken from the equa- 
tor towards each pole. Differences between the northern 
and southern hemisphere can be seen in the top panel. 
It is important to notice that local-helioseismology in- 
ferred flows at high latitudes have been shown to be af- 
fected by the periodic variation o f the solar Bn-an^le 
(iGonzalez Hernand ezet al. I 120061 : IZaatri et al~t 120061 : 
Beckers 2007). Although we do not yet have a full un- 
derstanding of the effect, it is suspected to be related to 
the foreshortening in the data associated with the vari- 
ation in the solar Bo-angle, the inclination of the solar 
rotation axis towards Earth. This periodic variation is 
clearly visible in the top panel of Figure [21 This system- 
atic effect definitely affects our results; however, it should 
be consistent throught the solar cycle and hence it can 
not be responsible for the observed amplitude variation. 
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Fig. 2. — Temporal variation of the fitted polynomial to the 
meridional circulation observations at a depth of 5.8 Mm. A sym- 
metrical plot averaging both hemispheres is shown in the bottom 
panel. Positive velocities are taken towards each respective pole. 

An asymmetry between the North and South hemisphere 
is present in the results, in particular the fact that at the 
shallow layers the flows deviate from zero at the equa- 
tor. The neutral line seems to be displaced towards the 
South during this period, which coincides with a more 
active period of the southern hemisphere (see Fig [3]). 

The bottom panel of Figure [2] shows a symmetrical plot 
of the flows, by averaging the North and South hemi- 
sphere. This reduces the systematic effect of the Bo- 
angle, without correcting it. The amplitude variation of 
the meridional flow, anticorrelated with the solar cycle, 
is more clearly visible here. 

The extra circulation, in the form of bumps, superim- 
posed on the smooth pattern of the meridional circula- 
tion, can also be seen in Figure [T] to a depth of ~ 11 Mm. 
For the earlier years, the bumps are centered in latitudes 
close to ~30°, towards the end of 2006, they have al- 
most disappeared, only to appear again at > 40° at the 
beginning of 2007. 

3. A MERIDIONAL COMPONENT OF THE TORSIONAL 
OSCILLATION? 

The residuals from a low order polynomial fitting to 
the inferred meridional circulation at a depth of 5.8 Mm 
are shown in Figure [H The top panel is a proxy to the 
corresponding surface activity, calculated using a longi- 
tudinal average of MDI synoptic magnetograms over the 
studied period. The central panel presents the residuals 
averaged over 60 days and using a bilinear interpolation 
(IDL, REBIN function) for imaging the results. The bot- 
tom panel show a symmetric version, averaging North 
and South hemispheres. Positive velocities in both cases 
are taken towards the poles. The variation of this extra 
component of the meridional flows with the solar cycle 
resembles, to some extent, the behavior of the so-called 
solar torsional oscillation of the zonal flows during this 
extended solar minimum (|Howe et al. Il2009[ ). It appears 
at medium latitudes 40-50°) around 2007 and seems 
to be moving towards lower latitudes with time. Is this 
a meridional component of the torsional oscillation? 

The appearance of the bumps at high latitudes before 
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Fig. 3. — Temporal variation of the meridional circulation resid- 
uals at a depth of 5.8 Mm (central panel). Positive velocities are 
directed towards the poles. A symmetrical plot averaging both 
hemispheres is shown in the bottom panel. The top panel shows 
the location and magnetic strength of the activity during the same 
period (calculated from MDI synoptic magnetograms). 

the magnetic activity surfaces happens at a similar time 
as the onset of the torsional oscillation, with the center 
between faster and slower bands in 2007 located around 
40° latitude. However, the amplitude of these residuals 
seem consistently smaller than those from the torsional 
oscillation, although marginally. Also, unlike the tor- 
sional oscillation, the extra component of the meridional 
circulation does not seem to converge towards the equa- 
tor at the end of the solar cycle. The bumps are also less 
clear below a certain depth and seem to disappear in the 
range attained with our analysis, although this maybe an 
artifact of the smoothing in the inversion technique. The 
results from global helioseismology analysis show that 
the torsional oscillation of the zonal flows remain present 
throughout the convection zone (Vorontsov et al. 2 QQ2[ 
IHowe et al. I |2005j); ring-diagram analysis shows little 
variation of the zonal flow amplitude in the outer 15 Mm. 
The center of mass of the perturbance seems to be stable 
above 20 degrees latitude. Even with our limited spatial 
resolution, if the bumps were located directly under the 
surface activity, the north and southern residual branches 
would be much closer to the equator towards the end of 
the cycle. 

Figure |4] represents an illustration of the meridional 
flow (top panel. A) with superimposed bumps, inwards 
flows towards the center of activity (top panel, B). We 
move the bumps with time following the latitude loca- 
tion of surface activity, keeping the amplitude of bumps 
constant, and then subsample the results with the same 
resolution that we obtain from the ring-diagram analysis. 
To image the results, we use the same bilinear interpo- 
lation function (IDL, REBIN). The central panel shows 
what happens when the activity is only present at the 
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surface. It can be seen that the residuals continue their 
convergence towards the equator after 2006, when they 
disappear from the real observations. To check if this 
was an effect of the new components at high latitudes, 
we repeat the experiment now including extra bumps 
centered at 50° in 2006 and moving slowly towards the 
equator. Neither of these two simple experiments agree 
with the observations. In the real data, the lower-latitude 
bumps are suppressed from 2006 to 2008, when surface 
activity is still present. However, a simple rel ation is 
not expected because, as found in previo us work (|Gizon I 
120031 : iGonzalez Hernandez et.~an I2008D . removino: the 
areas with surface activity changes the meridional cir- 
culation. Hence, there is also an expected contribution 
to the flows from active regions. Also, the coupling be- 
tween zonal and meridional flows through the Coriolis 
force is latitude dependent. 

Local helioseismology inferences of meridional circu- 
lation during the previous minimum do not show clear 
evidence of this medium-latitude formation of extra 
circulation before the activity appears at the s urface 
(|Chou and Dani200ll : [Beck. Gizon Duvah II2002I ). Un- 
fortunately, the availability of high resolution helioseis- 
mic data for local helioseismology methods dates back 
only to 1994, with the TON instrument, and to 1996 
with MDI. This, combined with the short span between 
the end of solar cycle 22 and the beginning of solar cycle 
23 and the long term averages needed to infer meridional 
circulation makes the isolation of the medium latitude 
branch difficult for that period. Thanks to this extended 
solar minimum we have been able to clearly observe the 
phenomenon before the onset of the new cycle. 

4. DISCUSSION AND CONCLUSIONS 

Our results confirm the solar cycle variation of the 
overall amplitude of the meridional circulation, minimum 
amplitude at maximum activity and vice versa. The 
trend is clearly there during the decline of solar cycle 
23 and has continued during the current prolonged solar 
minimum. No clear sign of trend reversal is found in the 
studied period. 

We also present, for the first time, a conclusive proof 
that the extra circulation, or bumps, of the meridional 
circulation under the solar surface is independent o f the 
presence of surface magnetic activity. 'Spruit ' ('2003') pre- 
sented a model that explains the torsional oscillation as a 
geostrophic flow due to the lower subsurface temperature 
in active regions. The model predicts the appearance of 
flows from the edges towards the center of the main lati- 
tude of magnetic activity, a meridional version of the tor- 
sional oscillation with a maximum amplitude of ~6m/s 
at the surface. The model also predicts a rapid decline of 
these oscillations with depth, which would disappear be- 
low 30 Mm. Our results are consistent in amplitude with 
this model. However, we show a more rapid attenuation 
of the inflows with depth, which disappear around 10- 
14 Mm and we are observing these flows at medium-high 
latitudes in the absence of surface activity. 

The extra component at medium-high latitudes starts 
to develop in 2007 for the new cycle, and continues there 
for the rest of the extended minimum. With our lim- 
ited resolution and data span, it is difficult to conclude 
if there is a variation in the latitude of the extra com- 
ponent for the last two years, although it seems to be 
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Fig. 4. — Illustration of the extra meridional circulation. Top 
panels show the assumed meridional circulation, the model bumps 
and the result when adding them both at about 20° latitude. The 
central panel presents the residuals, with the ring-diagram reso- 
lution, when moving the bumps acording to the location of sur- 
face activity. The bottom panel includes the developing branch at 
higher latitudes. These two simple experiments do not agree with 
the observations, at least close to the equator. 

migrating slowly towards lower latitudes. Applying ring- 
diagrams method to the current helioseismic data it is 
not possible to reach latitudes higher than ~50°, hence 
we can not distinguish the possible formation of this 
medium-latitude new branch earlier on at higher lati- 
tudes. The upcoming launch of the Solar Dynamics Ob- 
servatory (SDO) with the higher-resolution Helioseismic 
and Magnetic Imager (HMI) instrument on board will 
help to study meridional circulation above our present 
limits. 

The amplitude of the meridional residuals seems con- 
sistently smaller than those from the torsional oscilla- 
tion. A full, comprenhensive model that include both 
the zonal and the meridional components of the resid- 
ual circulation is necessary to understand the dynamic 
implication of this cycle varying component. 
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